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Abstract 



,-Q We study the survival probability of an immobile target in presence of N independent 

diffusing walkers. We address the problem of the Mean Target Lifetime and its depen- 
dence on the number and initial distribution of the walkers when the trapping is perfect 
or imperfect. We consider the diffusion on lattices and in the continuous space and we 
address the bulk limit corresponding to a density of diffusing particles and only one 
isolated trap. Also, we use intermittent motion for optimization of search strategies. 
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O , 1. Introduction 

,0, 

Several models of trapping of random walkers by a target have been extensively 
discussed in the literature [l|-|4|]. In most cases, in the analysis of the target lifetime, the 
^- ' initial distribution of walkers is taken at random (spatially homogeneous distribution) 

and the trapping is assumed "perfect" 5-7]. On the other hand, some results concerning 
to the order statistics of a set of random walkers are known when they are initially 
placed onto a given site of an Euclidean lattice [g, [9j or a fractal structure [2j. Thus, 
the question that arises naturally is how different initial configurations of the walkers, 
number of walkers, and trapping mechanisms affect on the survival probability of the 
^) ' target. 

In this paper we address the problem of trapping by a fixed target at the origin in 
the presence of a set of N independent random walkers. The quantity of our interest 
is the lifetime of the trap, which ends when any walker reaches the target under the 
appropriate circumstances. Our approach not only provides a unified framework that 
comprises several situations of trapping scattered in the literature, but also enables us to 
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compute (analytical and numerically) the lifetime of the trap in an exact way for a vast 
number of practical situations. Thus, we consider that the wandering of the walkers in the 
space may be normal or anomalous diffusion, and that the trapping may have particular 
characteristics: It may be perfect, in which case the lifetime of the trap reckons the time 
spend by the first walker to reach the origin, or imperfect [a, |lQj. The last situation 
includes the cases in which a walker passing by the origin is not trapped with certainty. 
Also, the present formalism could be applied even though we have dynamical trapping, 



i.e., the state of trap change in time like gated trapping [llL Il2j|. 

A novel application of our concepts comes from the hand of intermittent motion 
phenomena, which has recently motivated numerous studies in physics, chemistry, and 



biology |13l - ll7l |. In this piece of work, we show that the analysis of the trap's lifetime, as 
a function of the transition rates among internal states of the walker, allow us to optimize 
the intermittent search strategy for a hidden target. 

The paper is organized as follows. Section [5] presents the general formalism and 
define the concepts of survival of the target, lifetime density, and mean quantities; and 
establishes the connection with the problem of only one walker. Particularly, the last issue 
or first-passage time problem [18j . is reviewed in Sec. O whereas in Sec. [4] we reconsider 
the effects of dimensionality and number of walkers on our problem. This section also 
reviews the basic about the continuous-time random walk (CTRW). In Sec. [5] the effects 
due to initial spatial distribution of the set walkers are analyzed and the bulk limit is 
constructed. Section [5] presents several assorted illustrations for discrete and continuous 
systems with different types of initial distributions and considering the effects of finite 
size of the space and imperfection in the trapping mechanism. Last, in Sec. [3 we discuss 
the usefulness of our approach for searching targets with stochastic intermittent motion. 
Finally in SecJS]we give our conclusions. 

2. Lifetime of the Target 

The major objective of this contribution is the study of the effects of the initial 
distribution of independent walkers and the characteristics of the trapping process in the 
survival of the target. For this task, we begin reviewing and generalizing the formalism 
developed in Ref. [5|, |||. The survival probability at time t, $iv(i), of the static target 
(trap) at the origin in presence of N independent walkers that diffuse on a lattice can be 
written as [5j 

N 
s*i sn i= 1 

where u(s\, ..,sjv) denote the joint probability distribution of initially finding the first 
walker at a position s%, the second at sb and so on. $i(s*i, t) is the survival probability of 
the target at time t in the presence of only one walker initially at position Sj. In Eq. (TTJ), 
the sums run over all the lattices sites, and become in integrals over the space in the case 
of diffusion in the continuous space. 

From the survival probability, we define the target lifetime density (TLD), Fpi(i), in 
the presence of N walkers, in the standard way by 

F N (t) = ~* N {t). (2) 



Then, we can write [a, |6j 

N N 

F N (t) = J2 ■ ■ ■ Z)«(Si. ••' s ^) 2 Fl ^' *) II *i(^.*)i ( 3 ) 

si sjv i=l J#i 

where iq(sj,t) is the target lifetime density in the presence of only one walker, initially 
at position Sj. This quantity is defined from $i(s*j, i) in an analogous way to Eq. ([2)1. 

Fi(Si,t) = -^$i(Si,t). (4) 

In the case of perfect trapping, Fi(si,t) is the first-passage time density of the walker. 
When Fi(Si,t) is normalized, trapping is certain and the process is called recurrent in 
the sense proposed by Hughes [19(. On the other hand, if 



/i(si)= / F 1 (s ll t)dt<l, (5) 

Jo 

then the process is called transient [19| . fi(si) is the probability that a walker starting 
from site s*, will ever reach the origin. 

Now, using TLD, we introduce the Mean Target Lifetime (MTL) @,[20| 



t/v 



tF N {t)dt. (6) 



o 



If £ $j\r(4) —> for £ — > oo, then we can also write 



T N = / $ N {t)dt. (7) 

Jo 



3. First-Passage time 



A general expression for $i(si,£) can be constructed in terms of the conditional 
probability q(s,t\si,t = 0), corresponding to a walker be in s at time t, given that it was 
at Si at t — 0, restricted by the presence of a trap at the origin 

*i(Si,t) = 5^g(s;t|Si,t = 0) ) (8) 

s 

where the sum runs over all lattice sites and must be replaced by an integral in the 
continuous case. This expression is valid for any kind of trap, allowing for example 
imperfect trapping or dynamical gated trapping. 

For Markov processes, in the perfect trapping case, we can additionally exploit the 
connection between the probability density of first arrival at the origin at time t from the 
initial site s*i, -Fi (si,t), and the conditional probability of finding an unrestricted walker 
at site s at time t, given that it was initially at Sj, P(s, t\si, t = 0), [21| 

P{6,t\s ll t = 0) = tf(3i,t)<7 ? . g + / P(0,i|0,i')^iOM')^', (9) 



where ^(si,T) is the sojourn probability, i.e., the probability that the walker remains 
on the site Si a time lag r without a transition. Moreover, for an stationary process 
we have P(s, t\0, t') = P(s,t — t'\0,t = 0) and the integral in Eq. © becomes in a 
convolution. Thus, the Laplace transform of Eq. © lead us to the Laplace transform of 
the first-passage time density 

P(6,u\si,t = 0)-^(s i ,u)6 ffiiS 
l<i(si,u) — ^^ - , (10 

where the caret denotes the Laplace transform of the corresponding function. Therefore, 
using Eqs. (Ql and (TTU|) . the initial condition $i(si,i = 0) = 1, and taking Sj 7^ 0, we 
finally get the Laplace transform of the survival probability of the target in presence of 
only one walker, 

SxQM) = *fl*l g '* 7 ? -?&«&>* = °) . (ii) 

uP(0,u\0,t = 0) 

Finally, in a similar way as in Eq.Q, the mean first-passage time (MFPT) for the walker 
can be written as 

/>OG 

T = T 1 = / <f>i(t)dt. (12) 



4. MTL in d dimensions 



If the trapping is perfect, MTL is also the MFPT for the first of the set of N random 
walkers in reach the target. An interesting problem is presented when the MFPT for 
only one walker diverges, i.e., when the integral of Eq. (|12p diverges. In this situation, 
the interesting question that arises is to find the minimum number of walkers such that 
MTL becomes finite [8|, independently of the initial positions of the walkers. That is, 
from Eqs. ([T]) and (J7J, to find N for which the integral 

> N 

IJSifot)*, (13) 

t=i 

converges. Using our formalism, we can directly rederive the known results [8| for this 
problem. 

For concreteness, in this section, we use the CTRW [22|,|23| for the walker's dynamics. 
This allows us compute analytically the survival probability of the trap in presence of 
one walker. Hence, we can write 

P(s, u|0, t = 0) = ^ G(s, V-(u)) , (14) 

u 

where ij)(u) is the Laplace transform of the pausing time probability density, ip(t), and 
G(s, z) is the lattice Green's function. In d-dimensions, it is given by 



where A(k) is the structure function of the lattice. For a symmetrical walk on a simple 
cubic d-dimensional lattice we get A(k) = (cosfci + . . . + cos kd)/d, where fc, is the i-th 
component of k. 

Assuming that ip(u) — > 1 for u — > 0, the behavior of G(s, 4>{u)) is given by the values 
of k such that A(jfe) « 1, i.e., \k\ « 1. Thus, cos/c; « 1 - kf/2 and A(fc) = 1- |fc| 2 /(2cf). 
Therefore, we only need to consider the expansion exp(— is-k) sa 1 — is-k — (s-k) 2 /2 in the 
numerator of Eq. (|15p . Hence, for u — » 0, we get G(s,tp( u )) ~ G(0, -0(u)) — Kd(s,t/j(u)), 
where 

(s-k) 2 d d k ( 16 ) 



_* 1 - V(«) + ^(u) fc 2 /(2d) 

Given that we are dealing with unrestricted and spatially homogeneous walks, we get 
P(s,t\si,t = 0) = P(s — sl,t\0,t — 0). Hence, using these results in Eq. (jlip . we obtain 

$i(*,u)* * rf( *'f (u)) , (17) 

uG(0i(«)) 

where we have used Eq. (TH)) and the translational invariance of a CTRW on unbounded 
lattices (P(6,t\s,t = 0) = P(s,t|0,i = 0)). 

Kd{s, z) remains finite when 2=1 since the fc 2 in the numerator just cancels the 
singularity in the denominator. On the other hand, the singularity at z = 1 in G(0, z) 
depends on the dimensionality of the lattice. For studying this singularity, a good ap- 
proximate expression of Eq. (1151) is obtained taking the integrals over the d-ball of radius 
7r inscribed into the first Brillouin zone. Thus, for u — > 0, we get 

g (°^«))«7ct/ f-r f7 x l2 „ » » ( 18 ) 

(27r) a J l-ip( u )+ip(u)k 2 /(2d) 

where G d = 7r d / 2 /r(d/2 + 1). 

For normal diffusion, the pausing time density is ip{t) = Aexp(— At) [22j. Hence, 



Mu)=(l + j) , (19) 

and we get ip{u) ~ 1 — u/X, for u — ¥ 0. For anomalous diffusion, this asymptotic behavior 
is generalized by ip(u) s=y 1-K« a , with < a < 1. Hence, normal diffusion is the most 
severe case when we analyse the divergence of G(0, z) at z = 0. 

4-1- One dimension 

For d — 1, Ci — 2, and by direct integration of Eq. (TT51) we obtain 



G(g,^(u))w- 1 facbfmlirJ-^0- ) . (20) 



For u — > 0, ip(u)/{l — 4>(u)) — > oo and the last factor of Eq. (j20j) goes to tt/2. Thus, for 
normal diffusion, 

G(0,^(u))» 1 = = cxu~ 1/2 . (21) 

V2(l-te) 

In this manner, from Eq. (|17|l . we obtain $i(sj,u) ex it -1 ' 2 and using a Tauberian 
theorem 2fJ,|23|], we immediately get $i(s*i,t) c>c i -1 / 2 . Therefore, Eq. (fT2|) diverges and 



the convergence of Eq. (TT51) is obtained for N > 2. This result is also valid for imperfect 
traps [6J, [lOj] , and for dynamical gated trapping [llj, [l2| too. 



^.#. Two dimensions 

For d = 2, C2 = 7r, and by direct integration of Eq. (TT5|) we obtain 




G(0, &«)) « -^- In 1 + " ^ . (22) 



Thus, for m — >• 0, G(0,i/j(u)) ~ — ln(l — ip(u))/n ex — lnu, $i(s*j,u) ex — l/(u In it) and 
$i(si,i) ex hit. Therefore, the convergence of Eq. (|13[) is not reached for any value oi N. 

4.3. d > 3 

For d > 2, G(0,ip(u)) remains finite for u — > since the fc d_1 in the numerator just 
cancels the singularity in the denominator of Eq. (JT5J). Thus, from Eqs. (TT71) . $i(s*i, u) oc 
1/u. Alternatively this behavior can be seen from Eq. Q, and taking into account 
Eq. ([S]). Thus, the asymptotic behavior of the survival probability results 



$i(Si,t-> ex)) = 1 - / Fi{si,t')dt' = 1 - /!(£) 

Jo 



(23) 



Therefore, the long time behavior of the survival probability is constant, i.e., time inde- 
pendent. For d > 3, the process results transient [19] (0 < fi(si) < 1) so each factor in 
Eq. (fT3")) results < 1 — /l(fii) < 1, and then, the convergence of Eq. (J13D is not reached 
for any value of N. 

On the other hand, although the MFPT diverges for d > 3, the asymptotic limit 
of the survival probability of the target, which is proportional to IL_i(l — /i(^0) with 
each factor less than one, decreases monotonously with N , independently of the initial 
distribution of the walkers. 

Alternatively, we can proceed as in Ref. [8( choosing only the random walkers that 
will ever reach the target. For this purpose, it is necessary to define the conditional 
probability density of first arrival at the origin at time t from the initial site Si, given 
that the walker will eventually arrive there: Fi(si,t) / fi(si). Working with this quantity, 
we can find in d = 3 an expression equivalent to Eq. (1211) . Thus, in d — 3, MTL is finite 
if at least three walkers eventually reach the origin (target). 

5. Initial joint distribution 

Now, we consider the effects on TLD and MTL due to different initial probability 



distributions, u(si, .., s/v). Most of cases in the literature [24J426J belong to the following 
kinds of distributions: 

6 



5.1. Concentrated 

All walkers can begin at the same point of the space, 



N 



tt(Si,..,SV) = n*3i,«b, («o^O). (24) 

Thus, from Eq. ((TJ we obtain 



$rw = ($i(4t)) w , (25) 



and using Eq. © results 



F c N ° n (t) =NF 1 (? ,t) (^(Zo^f- 1 . (26) 

Recalling that $i(s*o,i) = 1 — L Fi(so,r)dr, the last equation may be compared with 
Eq. (19) inRef.Q. 

5.2. Equally likely sites 

Alternatively, the initial site of each walker can be chosen by chance among the M 
sites of a given set iS with equal probability, 

,- - x / M- N if all^eS, 
«(*W") = ( otherwise, (27) 

where S is such that ^ S. Thus, from Eq. (JTJ results 

JV 

(28) 



*#'(') = Iji E **(**) ) 



Also, we can write Eq. (|28p as &ff(t) = (($i(s, £))) , where (• • •) denotes the spacial 
average taken in the set S. Alternatively, we can recast Eq. (|2"5]l as 



*#*(*)= (i-^£(i-*i(*i,t))J 



(29) 



where Si (i — 1, . . . , M) are the positions of the M sites of set S. Note that there is 
not any restriction between N and M, Eq. (121)1) allows us to take the limits N — >• oo, 
M — > oo, with N/M — > /3 constant, i.e., the bulk limit. In this case we get 

<&f(i) = cxp(-/3S(t)), (30) 

where 

S(t) = 5^(1 -$i (?«,*)), (31) 

is the average number of lattice points visited, at least once time, by one walker, until 
time t. We are assuming that the series in Eq. (|3"Tj) converges. 

7 



In the continuous d-dimensional space, assuming that S has a finite volume V and 
that initially each walker begin uniformly distributed in V, the generalization of Eq. ((28|) 
is immediate, 

^ s (t)=(^J^i(s,t)d d s^ . (32) 

If the number of walkers per unit of volume, c = N/V, is constant, then in the limit 
TV — >• oo and V — > oo we also obtain 

$f {t)=exp(-c [(1 -*!(s, t))d d s) , (33) 

where we assume that the integral over the whole space in bounded. 

6. Illustrations 

6.1. Finite and semi-infinite chain 

As illustration, we now compute the MTL for a perfect trap at the origin of a chain, 
in presence of N walkers that jump from any site to its nearest neighbor with transition 
rate A. For a finite chain of L sites with absorbing end at the origin and reflecting end 
at site L, there is an exact expression for the Laplace transform of the first-passage time 



density [27 1 



R{ u y +R{u) 2L+] 
l + R{uf 



F i c?. «) = ; i nl \L+i — ' ( 34 ) 



(j = 1,2,..., L) where R(u) = (r + 1 — yjr 2 + 2 r) and r = u/X. This expression can 
be Laplace antitransformed in exact way [28j j but it is rather clumsy to display here. 
Moreover, for use Eqs. ([23)1 or (|28l) . we need previously make the integration involved 
in Eq. ([2]) and later make the integration in Eq. ([7]). We performed numerically these 
integrals [28j. In Fig. Q] we plot the values of T/v for a chain with L = 10 and initial 
distribution of walkers given by Eqs. (pM|) and (f2~T|). Strikingly, we find in both cases a 
power-law behaviour (see quasi-linear relation in the log-log plot) for almost all values 
of AT. 

If the chain is semi-infinite, an explicit and simple expression for $i(j, t) (j = 1,2,...) 
can be derive from results in the literature [y, l29j 

3 

*i(7,t) = e- xt (J (At) - J,- (At)) +2e- At ^J fe (At) , (35) 

fe=i 

where Ik(x) are the modified Bessel functions. Using this expression in Eqs. (|25|) or (j30|) 
and integrating numerically [28] in Eq. (J7J, we can evaluate the MTL for the initial 
distributions of Eqs. (|24l) and (|2T|) . Figure [5] plots the situation for walkers initially 
distributed with concentration /3 on the chain. 

For comparison purposes, we also include in Fig.[T]the plots corresponding to a semi- 
infinite chain with all walkers initially concentrated, at site s = 1 or site s — 10. In 
the later plots, the minimum number of walkers to obtain finite values of T/v is three, 
as has been quoted in Sec. 14.11 We want to stress the notorious resemblance obtained, 



for not so large values of N, between the cases of finite chain and semi-infinite chain for 
all walkers initially concentrated. Fig. [2] also includes the plots corresponding to a finite 
chain of L = 10 sites. In this graph the noticeable resemblance for large values of /3 is 
given between the cases of semi-infinite chain and finite chain with initial distribution 
of equally likely sites. 

6.2. Bulk limit in d- dimensions 

For d-dimensional lattices with a perfect trap at the origin, the bulk limit is given 
by Eq. ([50)1 . Using Eq. ([7]), MTL can be written in this case as 

exp(-/35(t)) dt. (36) 

It can be shown that the laplace transform of S(t), asumming a CTRW dynamics, is 
given by [30J 

S(u) = M , (37) 

u 2 P(0,u|0,i = 0) 

where ip( u ) ls given by Eq. (1191) . 

Using known expressions for P(0, u|0, t — 0) given in Ref. [19(, we plot in Fig. ((3]) Tp 
as function of (3 for different lattices in d = 1 (chain), d = 2 (honeycomb, square, and 
triangular), and d = 3 (SC, BCC, and FCC). Tp shows again a quasi-linear relation in 
the log-log plot and, as expected, Tp monotonously decrease as the coordination number 
of the lattice, k, is increased. 

6.3. Imperfect trapping in the continuous space 

The semi-straight line with an imperfect trap at the origin is another interesting illus- 
tration of our concepts. The probability density (in presence of the trap), q(x, t\x$, t = 0), 
for finding a particle at the location x at time t, given that it departed from site xq at 
the time t = 0, satisfies the classical diffusion equation 

where D is the diffusion coefficient. An imperfect trap is described by the radiation 

jq(x=0,t\x ,t = 0), (39) 



boundary condition 31 1 



dx 



x=0 



where 7 measure the efficiency of the trap. 7 = corresponds to a reflecting boundary 
and perfect trapping is reached in the limit 7 — >• 00. The solution of Eqs. (|3"5f and (|3"5|) 



can be read from Ref. 32] 



q(x,t\x ,t = 0) = q a -7%, (40) 

where 



(ar — Xq) \ / (x + xq 



cxp — + exp 



2Dt J "V 2Dt > 
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and 

q b = cxp (D-fH + j(x + x )) erfc ( ^±^£ + jy/Dt) . (42) 

V 2\ Dt J 

The relation between the survival probability $1(3:0, t) and q(x,t\xo,t — 0) is given 
by Eq. (J8]). Therefore, for the imperfect trapping we get [33[ 

*i(io, t) = 1 - erfc ( -^= J + exp (L> 7 2 t + 7x0) erfc ( ^= + -yVDt) . (43) 
\2y Dt J \2y Dt J 



Notice that the expected expression for perfect trapping (7 — > 00) is directly recovered 



Figure 0] graphs T/v for several situations of trapping efficiency (7) and for the initial 
distribution of walkers concentrated at Xq = 1. To make each plot we have used Eq. (|4"3"|) 
and (|2"5)l , and we have make numerically [28[ the integration of Eq. ([7]). Also in the 
continuous space, we obtain quasi-linear relations in the log-log plots of Tv vs N for 
not so large values of N. 

7. Target search with intermittent motion 

Examples of intermittent processes may be found in many fields. For instance, the 
case of a reactant that freely diffuses in a solvent and intermittently binds to a cylin- 
der |13| . Also, we found intermittent motion in the binding of a protein to specific sites 
on DNA for regulating transcription, as it is the case when the protein has the ability 
of diffuse in one dimension by sliding along the length of the DNA, in addition to their 



diffusion in bulk solution 14|. Moreover, intermitency could be associate with dynamical 
trapping problems [34|. We also found the search strategies like those implemented by 
animals in the pursuit of prey, or even in human activities such as victim localization, 
among the most representative examples of intermittent motion at macroscopic scales 
(see Ref. [15( and references therein). Recent works on intermittent search strategies 
were focused on the analysis of trapping of a single walker wandering in the presence 
of distributed traps. The magnitude that is usually estimated is this case is the named 
search time [171 ]. 

An interesting application of our present formalism is the related problem of a single 
static target among a set of initially uniformly distributed searchers, switching intermit- 
tently between two states of motion 15, 16]. In what follows, we calculate MTL for a 
target (perfect trap) at the origin of an infinite chain. We assume that each walker can 
be in either of two propagation states. In one of them, the displacement of the searchers 
is a random walk with symmetric jumps to nearest neighbors sites, with constant rate A. 
In the other state, the searchers also perform a symmetrical random walk, but jumping 
to next-nearest neighbors, with the same constant rate A. Thus, in the second state 
the diffusion is twice as big as in the other state, but in the second internal state, the 
walker can skip over the target without trapping. Transitions between the first and the 
second internal state, take place with rate constant 71, whereas the opposite transitions, 

10 



between the second and first state are at rate 72. The coupled master equations that 
describe this composite process of one walker are then 

^^ = \ (Pitt + l,t) + Pitt - M)) " XPi(j,t) +J 2 P2(j,t) - liPitt,*) , (45) 

^^ - \ (P 2 (j + 2, t) + P 2 (J - 2, t)) - A P 2 (j, t) + 7l Py(j, t) - l2 P 2 (j, t) , (46) 

where Pi(j,t) (P 2 tt,t)) is the joint probability that the walker be at site j with internal 
state 1 (2) at time t. If we assume equally likely sites as initial distribution of walkers 
and the bulk limit defined in subsection 15.21 we can use Eas. ((36l and (l37t with P(j = 
0, u) = Pi(j — 0, u) + P 2 (J = 0, it). In this manner, we can compute Tp as a function of 
71 and 72. 

In Fig. [5] we draw the behavior of MTL, in the bulk limit, as a function of the 
parameters of transition 71 and 72. As can be seen from the figure, we obtain a region 
of optimal values in the parameter space (71,72) which can be appreciated by the grey 
scale (darker means a smaller value in 1)3(71 ,72))- The valley in the surface indicates that 
we can tune the parameters to optimize the search. Notice how the MTL adequately 
characterize the improvement provided by intermittent search strategy [16j. Similar 
behavior has been exhibited in Ref. [17[ despite the difference in the addressed problems. 

8. Conclusions 

This paper provides a simple, general, and unified formalism for the lifetime statistics 
of a fixed target in presence of a set of independent hunters or gatherers that diffuse in 
the space. Our framework compresses normal and anomalous diffusion on lattices as well 
as in the continuous space. Also, our scheme allows us to consider perfect an imperfect 
trapping and can be directly extended for dynamical traps. Our main quantity, the 
MTL was introduced and its connections with any other physical quantities, relevant for 
the temporal statistics of trapping, were established. Particularly, the role of the initial 
spatial distribution was discussed. 

For trapping problems, where we deal with only one tramp, MFPT approach is limited 
when the process is transient (see Eq. ([5])) or the recurrence time in the lattice is infinite, 
because the MFPT diverge. Although MTL overcome this problem in one dimension, if 
the number of walkers is increased, it also diverges for d > 1. However, the bulk limit of 
MTL is finite in all situations. This robust property points out that MTL is the relevant 
physical quantity for situations where we consider the bulk density of diffusing particles 
and only one isolated capture center. 

Additionally, a striking feature of MTL is it presents non- universal scaling laws when 
it is plotted as a function of the number or density of walkers. Moreover, MTL is 
an efficient global optimizer for search strategies using intermittent motion. The MTL 
surface has a valley, allowing us to tune up the parameters that regulate the intermittency, 
and thus to minimize the time of search. 
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Figure and Captions 
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Figure 1: Tjv as a function of the number N of walkers for a chain with the trap at the origin. The 
finite case corresponds to L = 10 sites and two different initial distributions: Concentrated (CON), at 
site s = 1 or site s=10; and with all sites equally likely (ELS). The semi— infinite case (SI) corresponds 
to all walkers initially concentrated, at site s = 1 or site s = 10. The dotted lines are only to guide the 
eye. 



11 



100 - 



10- 



^T 



1 - 



. CL 


-V--. 




"~""^-v» ^tV""— V--- « 


-a tx ^ 








N A 




— D SI 




- 




^A^ 


- 


-o ELS 


~^A 


- 


-a CON,s=1 
-v-CON,s=10 


A ~A 


l t 1 1 1 1 1 1 1 — i — H 



1 2 3 45678910 

P(x10) 

Figure 2: Tp as a function of the concentration of walkers /3, for a chain with the target at the origin. 
The semi— infinite case (SI) corresponds to the bulk limit. The finite case corresponds to L = 10 sites 
(/? = N/L) and all walkers initially concentrated (CON), at site s = 1 or site s = 10; or with all sites 
equally likely (ELS). The dotted lines are only to guide the eye. 
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Figure 3: Tp as a function of the concentration of walkers, /3, in the bulk limit, for different lattices. The 
coordination numbers of the lattices, from top to botton, are k = 2, 3, 4, 5, 6, 8, 12. 
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Figure 4: Tjv as a function of the number JV of walkers for semi-straight line with an imperfect trap at 
the origin. The dotted lines arc only to guide the eye. 
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Figure 5: Tp in the bulk limit, for a concentration = 0.1 of walkers, as a function of parameters 71 
and 72. 
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